arXiv:1502.07759vl [astro-ph.SR] 26 Feb 2015 


Astronomy & Astrophysics manuscript no. CAbate © ESO 2015 

March 2, 2015 


Carbon-enhanced metal-poor stars: a window 
on AGB nucleosynthesis and binary evolution (I) 

Detailed analysis of 15 binary stars with known orbital periods. 

C. Abate 13 , O. R. Pols 1 , A. I. Karakas 2 , and R. G. Izzard 3 

1 Department of Astrophysics/IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands 

2 Research School of Astronomy & Astrophysics, Mount Stromlo Observatory, Weston Creek ACT 2611, Australia 

3 Argelander Institut fur Astronomie, Auf dem Hiigel 71, D-53121 Bonn, Germany 
e-mail: cabate@uni-bonn.de 

Received accepted ... 


ABSTRACT 

AGB stars are responsible for producing a variety of elements, including carbon, nitrogen, and the heavy elements produced in the 
slow neutron-capture process (^-elements). There are many uncertainties involved in modelling the evolution and nucleosynthesis of 
AGB stars, and this is especially the case at low metallicity, where most of the stars with high enough masses to enter the AGB have 
evolved to become white dwarfs and can no longer be observed. The stellar population in the Galactic halo is of low mass (< 0.85 M e ) 
and only a few observed stars have evolved beyond the first giant branch. However, we have evidence that low-metallicity AGB 
stars in binary systems have interacted with their low-mass secondary companions in the past. The aim of this work is to investigate 
AGB nucleosynthesis at low metallicity by studying the surface abundances of chemically peculiar very metal-poor stars of the halo 
observed in binary systems. To this end we select a sample of 15 carbon- and s-element-enhanced metal-poor (CEMP-i) halo stars 
that are found in binary systems with measured orbital periods. With our model of binary evolution and AGB nucleosynthesis, we 
determine the binary configuration that best reproduces, at the same time, the observed orbital period and surface abundances of 
each star of the sample. The observed periods provide tight constraints on our model of wind mass transfer in binary stars, while the 
comparison with the observed abundances tests our model of AGB nucleosynthesis. For most of the stars in our sample, we find that 
an episode of efficient wind mass transfer, combined with strong angular momentum loss, has occurred in the past. In some cases 
we find discrepancies between the observed and modelled abundances even if we adopt a fine-tuned set of parameters in our binary 
evolution model. These discrepancies are probably caused by missing physical ingredients in our models of AGB nucleosynthesis and 
they provide indications of how to improve our knowledge of the process of nucleosynthesis in AGB stars. 


1. Introduction 


In the final nuclear-burning stage of their lives, stars with initial 
masses between about 0.8 M 0 and 8 M Q ascend the asymptotic 
giant branch (AGB). During this phase of evolution the stellar ra¬ 
dius and luminosity increase by two to three orders of magnitude 
and the convective envelope of the star is expelled and enriches 
the interstellar medium with the products of stellar nucleosyn¬ 
thesis. AGB stars play an important role in our understanding of 
the origin of the elements (|Travaglio et al. 2004; Romano et al. 


2010. Kobaya shi et al.| |2011| >. Nuclear reactions in the interior 

of AGB stars are responsible for producing a variety of isotopes 
of elements, such as carbon, nitrogen, fluorine, sodium, magne¬ 
sium, and also elements heavier than iron produced by the slow 
neutron-capture process (Busso et al. 1999||He rwig, 2005 1 . 

Despite the importance of AGB stars for the chemical evo¬ 
lution of galaxies, several aspects of their evolution are not well 
understood. The physics of mixing is poorly constrained and the 
connected roles of convection, overshooting, rotation, and mag¬ 
netic fields need to be analysed in detail. Mass loss in AGB stars 
is highly uncertain, and most prescriptions in the literature are 
based on semi-empirical fits to observations of AGB stars in our 


Galaxy or in the Magellanic Clouds (e.g. Vassiliadis & Wood 
1993] |van Loon et al.| |2005[ >. Other uncertainties include nu¬ 
clear reaction rates and low-temperature opacities that follow the 


chemical composition of the AGB model in detail; both are cru¬ 
cial for determining the predicted level of chemical enrichment 
from an AGB model (e.g.|Izzard et al.|[2007; Marigo & Aringer 
|2009] >. 


AGB nucleosynthesis depends on all these physical quan¬ 
tities, and therefore the observation of chemical abundances in 
AGB stars and their progeny (e.g. post-AGB stars and planetary 
nebulae) are an important source of information on the physi¬ 
cal processes that drive the evolution of these objects (|Herwig| 


2005). However, because the duration of the AGB phase is typ 
ically less than 1% of the total stellar lifetime (Vassiliadis & 


|Wood[|T993| , AGB stars are relatively rare compared to stars in 
earlier evolutionary stages. This is a limit especially at the low 
metallicity of the Galactic halo, where the stellar population is 
about ten billion years old and where AGB stars more massive 
than approximately 0.85 M G have already become white dwarfs. 


One way to overcome this limit is to look for signatures of 
AGB nucleosynthesis in halo binary stars. Consider a binary sys¬ 
tem where, in the long distant past, the primary star ascended the 
AGB. When this happened, some fraction of its stellar wind was 
transferred onto the companion, a lower-mass secondary star. If 
the total mass of the secondary star after the accretion did not 
exceed approximately 0.85 M Q , this star has not yet become a 
white dwarf and can in principle still be observed today. 
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This binary scenario has been invoked to explain the peculiar 
chemical abundances determined in the carbon-enhanced metal- 
poor (CEMP) stars widely observed among the very metal-poor 
(here defined a'-Hj [Fe/H] < -2.0) stars of the Galactic halo. 
CEMP stars constitute a significant fraction of metal-poor stars 


in the halo, between 9% and 25% (e.g. Marsteller et al. 

2005 

Frebel et al. 2006[ |Lucatello et al. |2006[ jCarollo et al. 

2012 

Lee et al. 2013), with the CEMP frequency rising witl 

a de- 


et al. 2013|). Most CEMP stars are also enriched in heavy el 


ements produced by slow and rapid neutron-capture processes 
(s-process and r-process, respectively). Traditionally the ex¬ 
cess of barium in stars is related to .v-process and the excess 


of europium to r-process (Sneden et al. 20081. Jonsell et al. 
( 2006| ) indicate as CEMP those very metal-poor stars with ob¬ 
served [C/Fe] > 1 and as CEMP-.? the CEMP stars that sat¬ 
isfy the relations [Ba/Fe] > 1 and [Ba/Eu] > 0 simultaneously. 
CEMP-.s stars with [Eu/Fe] > 1 are defined C’EMP-.s'/r stars. 
Other authors adopt slightly different definitions (e.g.,|Beers & 


|Christliebl|2005[|Aoki et afT] |2007[ |Masseron et al.| 201 0| > . The 

binary formation scenario is supported by detection of radial ve¬ 
locity variations in a large number of CEMP-? stars, statistically 
consistent with the hypothesis that all CEMP-s stars are in bi¬ 
nary systems ( Lucatello et al.[ 2005; Starkenburg et al. |2014| i. 

Several authors have used AGB nucleosynthesis models with 
the aim to reproduce the abundances observed in CEMP-.? stars 


(|Stancliffe & Glebbeek 2008 Stand iffe] 2009 [ |Bisterzo et al. 


2009 Masseron et al. 26l0||Bisterzo et ak 2011 2012 Lugaro 


2012||Placco et al.[ |2013| ). In these studies the observed 


et al. 


abundances of the CEMP-? stars are directly compared to the 
outcome of detailed models of AGB nucleosynthesis of vari¬ 
ous masses and metallicities, with some assumptions to estimate 
the dilution of the accreted material in the envelope of the sec¬ 
ondary star. In particular, Bisterzo et al. ( 2012| > provide an in¬ 
dividual analysis of 94 C 
-3.6 < [Fe/H] < 


MP-? stars in the metallicity range 
1.0. Placco et al. ( 2013| > have recently 
performed a similar analysis of two newly discovered CEMP-s 
stars for which they had previously determined the surface abun¬ 
dances of 34 elements. Bisterzo et al. (2012) and Pl acco et al.| 
(2013) conclude that the observed abundances of most CEMP- 
.s stars are consistent with the hypothesis that these stars have 
accreted mass from an AGB companion, while the models of 
AGB nucleosynthesis do not agree well with the abundances of 
CEMP-s/r stars. The physical process that leads to simultaneous 
enhancement of s- and r- elements in these stars is still debated. 


In this paper we focus our analysis on a sample of 15 ob¬ 
served CEMP-.s binary stars with known orbital periods. We 
tackle the problem of their formation history by studying the 
whole dynamical and chemical evolution of the binary systems. 
The purpose of this work is to study under which conditions our 
model reproduces at the same time the evolutionary stage, chem¬ 
ical abundances and orbital period of each observed star, and 
to understand which constraints can be placed on our model of 
binary evolution and of the nucleosynthesis processes in AGB 
stars. The measurement of the orbital period, besides implying 
that these systems have very likely interacted in the past, pro¬ 
vides a constraint on the initial orbital separation and hence also 
on the initial masses of the stars in the binary system. A larger 
sample of CEMP-.? stars without measured period is analysed in 


1 [X/Y]= log 10 (We/Ay) - log 10 (A X o/X Yo ), where A x ,y indicates the 

number density of the elements X and Y, and O denotes the abundances 
in the Sun. 


our forthcoming paper (Abate et al. 2015 submitted, hereinafter 
Paper II). 

The paper is structured as follows. In Sect. [2] we describe our 
sample of observed CEMP-.s stars. In Sect. |3]we briefly discuss 
the evolution of a star in the AGB phase according to detailed 
models. In Sect. Q] we summarise the main characteristics of our 
binary evolution model, we explain how we included therein the 
results of detailed models of AGB nucleosynthesis and we de¬ 
scribe our method to find the best fit of the observed abundances 
of each star in our sample. Sect. [5]is dedicated to the comparison 
between the data and the outcome of our model. In Sect. [6] we 
discuss the results while Sect.[7]concludes. 


2. Data sample 

To perform our analysis we select a sample of 11 binary stars 
according to the following criteria: (i) measured orbital period; 
(if) iron abundance -2.8 < [Fe/H] < -1.8; (Hi) enhanced 
carbon and barium abundance, respectively [C/Fe] > 1 and 
[Ba/Fe] > 0.5. The restriction on the abundance of iron is mo¬ 
tivated because our model of AGB nucleosynthesis is tailored to 
reproduce the abundances at metallicity Z = 10 4 , roughly cor¬ 
responding to [Fe/H] ~ -2.3 (see details in Sects. [3] and [4 1 . We 
ignore systems in which only upper or lower limits are available. 
We add to this sample four more systems: CS22956 - 028, in 
which barium is only weakly enhanced, [Ba/Fe] = 0.38, but the 
strontium abundance is [Sr/Fe] = 1.39; CS29497-034, which is 
enhanced in carbon and barium ([C/Fe] = 2.69, [Ba/Fe] = 2.12) 
and has many observed elements but the iron abundance is low, 
[Fe/H] = -2.96; HD198269 and HD201626, in which barium 
has not been measured but the abundances relative to iron of 
other .?-elements, e.g. lanthanum, cerium and lead, are enhanced 
by more than 1 dex. In Table [T] we list for all 15 stars in our 
sample the observed orbital period, P or b, eccentricity, e, mass 
function, f(M c ) surface gravity, effective temperature, 7/o', iron 
abundance and number of elements observed. For every star we 
collect from the literature the absolute abundance of each ele¬ 
ment, 

A X = 12 + log 10 ^ , (1) 

where A x and Ah are the number densities of element X and hy¬ 
drogen, respectively. To compute the abundances relative to iron, 
[X/Fe], we use the solar abundances as determined by |Asplund| 
|etak1 ( |2009l ). 

In most of the stars of our sample the abundances determined 
by different authors are consistent within the observational un¬ 
certainties. In these cases we adopt the arithmetic mean of the 
absolute abundances and the maximum published uncertainty. 
In star CS22948 - 027 there are large discrepancies between the 
abundances published by different authors, up to 1 dex or more. 
These discrepancies are mostly due to the different atmospheric 
paramet ers adopted: [Preston & Sneden ( |2001| ) and |Aoki et ST] 
(2002b find a low surface gravity, log 10 (g/cm s ' 1 2 ) = 0.8 and 1.0 
respectively, wher eas Hill et al. ([2000), Barbuy et al. ( 2005| ) and 
Aoki et akf (2007 1 find higher values, log 10 (g/cm s -z ) = 1.8, 1.8 


and 1.9, respectrvely. The most recent abundances published by 
Barbuy et al. (|2005|) and Aoki et al. (|2007[i are obtained from 


spectra at higher resolution and signal-to-noise ratio compared 
to previous studies and include corrections for effects due to non¬ 
local thermodynamic equilibrium, thus in our study we adopt the 
average of their abundances. In stars HD198269 and HD201626 
the iron abundances published by Vanture| ( 1992b) are 0.8 dex 
higher than the values indicated by Van Eck et al. (2003), al¬ 
though the abundances of most ^-elements are consistent within 
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Table 1 . Summary of the orbital parameters, surface gravities, temperatures and chemical properties observed in our sample CEMP- 
s stars. 


ID 

P orb/d 

e 

f(M c )/M e 

logiofe/cms z ) 

Teff/K 

elem. 

[Fe/H] 

[C/Fe] 

[Ba/Fe] 

Reference 

BD+04°2466 

4593 

0.286 

0.076 

1.8 ±0.2 

5032 

20 

-2.1 

1.3 

1.6 

19, 26, 27 

CS22942-019 

2800 

0.1 

0.036 

2.2 ±0.4 

4967 

18 

-2.7 

2.2 

1.8 

6, 8, 10 

CS22948-027 

426.5 

0.02 

0.004 

1.8 ±0.4 

4800 

21 

-2.5 

2.2 

2.0 

17, 22 

CS22956-028 

1290 

0.22 

0.076 

3.9 ±0.1 

6900 

11 

-2.1 

1.9 

0.4 

5, 13 

CS22964-161A 

252.48 

0.66 

0.147 

3.7 ±0.2 

6050 

21 

-2.4 

1.6 

1.4 

25 

CS22964-161B 

252.48 

0.66 

0.216 

4.1 ±0.4 

5850 

22 

-2.4 

1.4 

1.3 

25 

CS29497-030 

342 

0 

0.002 

4.0 ±0.5 

6966 

33 

-2.5 

2.4 

2.3 

5, 13, 16, 18, 23 

CS29497-034 

4130 

0.02 

0.060 

1.65 ±0.15 

4850 

20 

-3.0 

2.7 

2.1 

3, 17, 21 

CS29509-027 

196 

0.15 

0.001 

4.2 ±0.1 

7050 

5 

-2.1 

1.5 

1.3 

4, 13 

HD 198269 

1295 

0.09 

0.107 

1.30 ±0.25 

4800 

9 

-2.2 

1.7 

- 

1, 14, 20 

HD201626 

407 

0 

0.075 

2.25 ± 0.25 

5200 

10 

-2.1 

2.1 

- 

1, 14, 20 

HD224959 

1273 

0.179 

0.096 

1.95 ±0.25 

5050 

15 

-2.1 

1.8 

2.2 

1, 14. 20, 28 

HE0024-2523 

3.14 

0 

0.049 

4.3 ±0.1 

6625 

17 

-2.7 

2.1 

1.6 

10, 11, 12, 15 

HE0507-1430 

446 

- 

- 

0.8 ±0.1 

4600 

5 

-2.4 

2.8 

1.3 

24, 29 

LP625^14 

> 4383 

- 

- 

2.65 ± 0.3 

5500 

31 

-2.8 

2.3 

2.8 

2, 3, 7, 8,21 


References. (U Vanture (1992b 

cl; (2) 

Norris etaTj( 1997); (3) Aoki et al. (2000); (4) Hill et al. (2000); (5) Preston & Sneden (2000); (6) Preston 

& Sneden(200 

:>; (7) 

Aoki et al. 

(2001 

; (8) Aoki et al. 

(20023; (9) Aoki et al. (2002c); (10) Cohen et al. (2002); (11) Carretta et al. (2002); (12) 

Lucatello et al. 

(2003 

); (13) Sneden et al. (2003); (14) 

Van Eck et al. 

2003); (15) Cohen et al. (2004); (16) Sivarani et al. (2004); (17) Barbuy 

CD 

O 

O 

LA 

8) Ivans et al. (2005); (19) Jorissen et al. ( 

2005); (20) 

Lucatello et al. 

(2005); (21) Aoki et al. (2006); (22) Aoki et al. (2007); 

(23)[Johnson et al. (2007); (24) Beers et al. (2007); (25) [T 

lompson et al. (2008); (26) 

Pereira & Drake (2009); (27) Ishigaki et al. <2010fT(28) 


the observational uncertainties. In the analysis of these two stars 
we use the abundances of iron and ,v-elements published by Van 
Eck et al.] (20031 that are obtained from spectra at higher reso¬ 


lution. Because the abundances of carbon, nitrogen and oxygen 
are not provided by |Van Eck et al.| ( |2003| we adopt the values 
indicated by Vanture ~{ 1992a|b) >. 

The minimum uncertainty that we assume is 0.1 dex in the 
chemical abundances and log g. For the effective temperature we 
adopt an uncertainty of 100K unless differently stated. Orbital 
periods and eccentricities of the 15 binary stars in our sample are 
mostly published without errors; in Sect. |4.3| we explain how we 
deal with this missing information and which parameters we ac¬ 
tually use in our study. For stars HE0507 - 1430 and LP625 - 44 
no information about the eccentricity or mass function is avail¬ 
able. We adopt / J orh = 12 years as period of LP625-44, although 
this is in fact a lower limit (|Aoki et al. 2002a)> and the real period 
is probably much longer (Hansen et al,||20f2) . 


3. Nucleosynthesis during the AGB phase 

For a detailed review on AGB evolution we refer to Herwig 
( 2005| >. Here we summarise the thermal-pulse cycles that take 
place in AGB stars. A thermal pulse (TP) occurs when the he¬ 
lium shell ignites and expands the outer layers extinguishing the 
hydrogen shell. Third dredge-up (TDU) may occur after the TP, 
when the the inner edge of the convective envelope moves in¬ 
ward (in mass) and mixes to the surface the products of inter¬ 
nal nucleosynthesis. During the interpulse period hydrogen is 
burned quiescently and the ,v-process occurs in the intershell re¬ 
gion where helium is abundant and (a, n ) reactions can be effi¬ 
ciently activated to produce free neutrons. 

In low-mass stars up to 3M 0 the main neutron source is the 
l3 C(a, «) 16 0 reaction. I3 C is formed in the intershell region by 
partial mixing of protons from the envelope at the end of a TDU 
event. The protons are captured by the 12 C to form a layer rich 
in 13 C, the “ I3 C pocket”. During the subsequent interpulse pe¬ 
riod the 13 C(ar, n) I6 0 reaction takes place in the 13 C pocket and 
releases free neutrons (Straniero et al.|[T995|l. One of the largest 


uncertainties in the study of the ,v-process nucleosynthesis is re¬ 
lated to the numerical treatment of the 13 C pocket: its mass is 
essentially a free parameter in AGB models but some constraints 
can be put by comparison to the chemical composition observed 
in planetary nebulae, post-AGB stars and CEMP-s stars. We re¬ 
fer to Busso et al. (1999]) for a thorough discussion of the nucle¬ 
osynthesis in AGB stars. 

In the models of |Karakas (2010 K10 hereinafter), computed 
with the Monash/Mount Stromlo (hereinafter Stromlo) code, 
the evolution of the stellar structure is computed as a first step 
and subsequently detailed nucleosynthesis calculations are per¬ 
formed using a post-processing algorithm (Lugaro et al. 2004; 
Karakas & Lattanzio| 2007|l. In the post-processing algorithm, 


protons are mixed into the intershell region by artificially adding 
a partial mixing zone (PMZ) at the deepest extent of each TDU. 
This method is described by Lugaro et al. (2004 and is similar to 
that of Goriely & Mowlavi ( 2000[ ). Lugaro et al. ( 201 2) explore 
the effect of different PMZ masses in the range [0, 0.004] M 0 
on the nucleosynthesis of AGB stars of different initial mass. 
The mass of the PMZ is non-zero only in low-mass AGB mod¬ 
els (M < 3 M 0 ; although as a test the case M = 5.5 M Q and 
Mpmz = 5 x 10 4 AL.) is considered). At higher mass (and 
Z < 10 4 ) the ingestion of protons into the helium-flash-induced 
convection zone combined with the high temperature at the base 
of the convective envelope leads to large energy production. This 
may significantly affect the structure of the star and such an ef¬ 
fect cannot be taken into account in the post-processing algo¬ 
rithm (Sect. 2.2 of K10). 


|Lugaro et al.| ( |2012| ) distinguish four regimes of neutron- 
capture process in their models of low-metallicity (Z = 10 4 ) 
AGB stars and each regime dominates in a different mass range. 
In models of mass above 3 M 0 the 22 Ne neutron source domi¬ 
nates and lighter v-elements (e.g. strontium) are favoured with 
respect to heavier ,v-elements (e.g. barium and lead). The 13 C 
neutron source dominates in models of mass below 3 M a . For 
masses between 1.75 M 0 and 3 M 0 13 C burns completely in ra¬ 
diative conditions, neutrons are captured in the thin l3 C layer 
and consequently the heavier .^-elements are favoured because of 


3 
































































































C. Abate et al.: CEMP stars: a window on AGB nucleosynthesis and binary evolution (I) 


the large number of neutrons per iron seed. On the contrary, for 
lower masses 13 C burns convectively, neutrons are released over 
the whole intershell region and thus the neutron-to-iron ratio is 
lower compared to the case when 13 C burns radiatively. In mod¬ 
els with masses below 2.5 M 0 , ingestion of protons associated 
with the first few thermal pulses produce 13 C nuclei, which burn 
convectively and .v-elements are produced even in the models 
without PMZ. This regime dominates for masses below 1.5 M G . 

Lugaro et al. ( j2012 1 calculate detailed nucleosynthesis over 
a grid of 16 initial masses between 0.9 and 6M G for Z = 10 4 . 
These include abundances predictions for 320 isotopes from hy¬ 
drogen through to 21() Po as a function of interior mass and time 


for each single stellar model. In Sect. 4.2 we describe the method 
we used to implement these results in our model of binary pop¬ 
ulation synthesis. 

4. Model of binary evolution and nucleosynthesis 

To study the chemical compositions of CEMP-.? stars, we 
use the binary evolution and population synthesis code 
binary_c/nucsyn described by|Izzard et al.| ( 2004| 2006[> and 
recently updated by Izzard et al.| ( [2009 201 0| > and Abate et al. 
( |201 3| >. This code combines a binary-evolution model based on 
the rapid binary stellar evolution prescriptions of Hurley et al. 
(|2002|> with a synthetic nucleosynthesis model developed by 
|Izzard et aL] ( |2004||2006||2009| l. In the first part of this section we 
briefly describe the most important characteristics of our model 
(Sect. |4. 1 [ >, In the second part we explain how we updated our 
synthetic nucleosynthesis model with the results of our most re¬ 
cent detailed models of AGB nucleosynthesis (Sect. |4~2] >. In the 
last part we describe the method used to determine the best fit of 
the abundances observed in the stars of our sample (Sect. [43]). 


Eq. (6) by (Boffin & Jorissenj ( |T988j ) with ctehl = 10, 
whera^lafBHL is a numerical constant normally between 1 
and 2. 

The purpose of comparing these two model sets is to under¬ 
stand if our models are limited in reproducing the observa¬ 
tions by our treatment of the mass-transfer process. 
Thermohaline mixing is assumed to be efficient: the ac¬ 
creted material mixes instantaneously with the stellar enve¬ 
lope. Stancliffe et al. (|2007j) suggest that this approximation 


is typically reasonable. In some cases we relax this assump¬ 
tion and simulate the conditions of highly inefficient thermo¬ 
haline mixing, in which the accreted material remains on the 
stellar surface until mixed in by convection. 
Common-envelope evolution is computed according to the 
prescription of Hurley et al. ( |2002 1 with a free parameter for 
the common-envelope efficiency of ejection set by default to 
acE = 1- We do not include accretion during the common- 
envelope phase. 

As initial composition of isotopes up to 76 Ge we adopt the 
abundances predicted by the chemical evolution models of 
Kobayashi et al. (j2011)> for solar neighbourhood stars at 


[Fe/H] ~ -2.3. For heavier isotopes that have not been cal¬ 
culated by Kobayashi et al. ( 201l) >, we adopt the solar distri¬ 
bution of abundances by Asplund et al. ( 2009} > scaled down 
to metallicity Z = 10~ 4 . 

The nucleosynthesis algorithms that compute the evolu¬ 
tion of the surface abundances through the first and second 
dredge-ups are based on the work by Karakas et al. ( 2002j > 
and Karakas & Lattanzio (2007!. A prescription for hot- 


bottom burning is also included for star more massive than 
approximately 2.75 M G . For more details we refer to Izzard 
|et al.| ( [2004l|2006||2009| ). In the next section we describe the 
treatment of the third dredge-up. 


4.1. Input physics 

We describe here the most important input parameters that need 
to be set in our model of binary evolution and nucleosynthesis. In 
this section we list all the options and in Sect.[5]we specify which 
assumptions are made in our analysis of each binary system. 


The wind mass-loss rate up to the AGB phase is parame- 
terised according to the Reimers ( |1975| ) formula multiplied 
by a factor // = 0.5 on the first giant branch. During the AGB 
phase the formula of Vassiliadis & Wood ( |1993[ > is used, 
with minimum and maximum values of the wind velocity 
v' w = 5 km s and 15l<ms respectively. 

The variations of angular momentum because of mass loss 
and the efficiency of the wind mass-transfer process are cal¬ 
culated with two alternative model sets, namely: 

• Model set A (default): the angular momentum carried 
away by the expelled material is computed assuming a 
spherically symmetric wind (Eq. 4 of | Abate et al.[|2013| . 
The wind accretion efficiency is calculated according to 
a model for wind Roche-lobe overflow (WRLOF) that 
includes the dependence on the mass ratio of the binary 
system, as in Eq. (9) of Abate et al. (|2013 i 


Model set B: a model of efficient angular momentum loss 
is adopted in which the material lost from the binary sys¬ 
tem carries away a multiple y = 2 of the average spe¬ 
cific orbital angular momentum (Eq. 2 of llzzard et al. 


|2010| and Eq. 10 of |Abate et al.|[20l3] ). An extremely ef¬ 
ficient wind mass-transfer process is simulated by adopt¬ 
ing an enhanced version of the canonical Bondi-Hoyle- 
Lyttleton prescription for the wind accretion rate, namely 


4.2. Third dredge-up nucleosynthesis 

The efficiency of third dredge-up is defined by: 

, _ M dv 
~ AM h ’ 


( 2 ) 


where Mpp is the mass dredged up from the intershell region 
and AMh is the core-mass growth due to hydrogen burning dur¬ 
ing the previous interpulse period. Thus, over a whole interpulse 
period the core grows by: 


AM c = AM h - M DU = (1 - d)AM H . (3) 


A TDU episode occurs with efficiency A when M c exceeds a 
threshold mass M c m ; n . The values of A and M c min are functions 


of mass and metallicity fitted to the detailed models of Karakas 


et al. 

(2002), Karakas & Lattanzio 

2007) and K10, according to 

the algorithms explained by Izzarc 

et al. (2004 see their Sect. 

3.4 for details). Izzard et al. (2009 

introduced a free parameter 


in the model, the minimum envelope mass for TDU M env m i n , to 
study the effects of efficient TDU at masses down to 0.8 M 0 . To 
better reproduce A as a function of time from the detailed models 
of K10 we set. 


Me, 


M, 


o 


( 0.15 

= i 0.88 M, - 0.73 
0.37 


if M, < 1 M q , 

if 1 M g < M* < 1.25 M q , 

if M, > 1.25 M q . 


2 In the paper by Boffin & Jorissen ( 1988 ' q- B hl this constant is indi¬ 
cated simply as a 
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M, =0.9 M q M* =1 M 0 M, =1.5 M @ 



Thermal pulse number 



Fig. 1 . Dredge-up efficiency A (top) and the abundances of carbon and barium (middle and bottom, respectively) computed at each 
thermal pulse with the binary_c and Stromlo codes (solid and dotted lines, respectively) for single star of mass (from left to right) 
M» = 0.9, 1, 1.5 M q and M PUZ = 2x 10 ~ 3 M 0 . 


With these prescriptions, single stars of initial mass above 
0.9 M q experience some TDU at metallicity Z = 10 4 , in ac¬ 
cordance with the models of K10. 

Each TDU modifies the surface abundances of the model 
AGB star, because material from the stellar interior that has been 
subject to nuclear processing is mixed into the envelope. For this 
reason it is essential to know the chemical composition in the in¬ 
tershell region at every thermal pulse if we are to reproduce sur¬ 
face composition of the detailed model with our nucloesynthesis 
algorithm. In our model the abundances in the intershell region 
of an AGB star of metallicity Z = 10 4 are stored in a table as a 
function of three quantities: the mass of the star at the first ther¬ 
mal pulse, the thermal pulse number and a free parameter, Mpmz, 
that describes the dependence of the chemical composition on 
the mass of the partial mixing zone according to the Stromlo de¬ 
tailed models. Our table includes all the isotopes provided by 
the models of Lugaro et al. ( |2012| . During the evolution of an 
AGB star of mass M„, at each TDU an amount of mass Mdu of 
the intershell region is instantaneously mixed in the convective 
envelope. The chemical composition of the dredged-up material 
is calculated by interpolating values of stellar mass and Mpmz in 
our table. As an example of the chemical evolution of the surface 
of an AGB star due to the TDU episodes, in panels d-i of Fig. 
[I] we show the abundances of carbon and barium as computed 
at each thermal pulse with the binary_c and Stromlo codes for 
stars of initial masses 0.9, 1 and 1.5 M 0 with M pmz = 2 x 10 5 . 

In Fig.[TJi we show the value of A as a function of the thermal- 
pulse number computed with our model and with the Stromlo de¬ 
tailed code for a single star of mass 0.9 M e . The detailed model 


predicts no TDU between pulses number 4 and 7 and some TDU 
up to pulse 19 with a maximum A ~ 0.14. In our model we have 
increasingly efficient TDU episodes up to pulse 13, when the 
envelope mass becomes lower than M envjnln = 0.15. Despite the 
difference in /I the total amount of mass dredged up in the two 
models is similar (0.012 M 0 and 0.014 M 0 with the Stromlo and 
binary_c codes, respectively), as shown in the bottom panel of 
Fig. 0 Consequently, also the surface abundances predicted by 
two models are similar (Figs. |T]:1 and [TJr ). Fig.[TJ) shows the value 
of /I for a single star of mass 1 M 0 . The Stromlo code predicts 
only two TDUs, a weak one after pulse number 5 and a stronger 
one after pulse number 6. This peculiar evolution is only ob¬ 
served in the model of a 1 M 0 star, but not for higher masses (see 
K10 for details), and is probably related to the dependence of 
the TDU phenomenon on the numerical treatment of convective 
boundaries, as discussed e.g. by Frost & Fattanzio| ( |1996| l and 
Mowlavi (1999). At mass M* = 1 M 0 we do not force our model 
to reproduce A of the detailed model and as a consequence the 
total mass dredged up with our model is almost 10 times larger 
(0.02 Mq rather than 0.002 M 0 ). Fig. [TJ: is the same as Figs.|lji 
and |T|i but for a star of mass 1.5 M 0 . In this case the main dis¬ 
crepancy is due to the fact that the detailed model predicts more 
TDU episodes and therefore the total mass dredged up in our 
model is smaller (bottom panel of Fig. |2j. However, the final sur¬ 
face abundances of carbon and barium are essentially the same 
(as shown in Figs. |T]f and[lj). At higher masses the two codes are 
in good agreement, as shown in Fig. [2] where we plot the maxi¬ 
mum value of the TDU efficiency, d max , and the total amount of 


5 


























































C. Abate et al.: CEMP stars: a window on AGB nucleosynthesis and binary evolution (I) 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 


Initial mass / M ; 

Fig. 2. Maximum TDU efficiency /l max (top) and totai mass 
dredged up M DUj tot (bottom) for stellar masses in the range 
[0.9,6.0] as computed with the binary_c (solid line) and 
Stromlo codes (dotted line). 


choose Ami = 34, Am 2 = 28, N a = 30, Apmz = 10. We consider 
circular orbits in all our models. 

The initial parameters are chosen as follows: 

• M\ varies in the range [0.9,6.0] M 0 . Up to 3 M 0 the grid 
spacing is AM\ =0.1 and AMi = 0.25 otherwise. 

• AG is uniformly spaced between 0.2 M 0 and 0.9 M 0 (AMi = 
0.025 .V/,). 

• a varies in the range [10 2 ,10 5 ] R 0 . The distribution of sep¬ 
arations is uniform in log 10 a (Alog io a/R 0 = 0.1). In the 
mass range considered here, stars at wider separation do not 
interact in our models. All stars in our grid are formed in 
binary systems. 

• Mpmz is always zero in AGB stars of mass M > 3 M 0 , other¬ 
wise we use the following values: 0, 10' 4 , 2 x 10' 4 , 5 x 10' 4 , 
6.66xl0' 4 , lxlO' 3 , 1.5xl0' 3 ,2xl0~ 3 ,3xl0'\ 4xl0' 3 M 0 . 


We evolve the systems on this grid in parameter space and 
with different combinations of the physical parameters described 
in Sect. 4.1 We select the stars that have not yet become white 
dwarfs at the evolutionary time t > 10 Gyr and therefore satisfy 
the surface gravity condition log 10 (g/cm s' 2 ) < 5.0. We compare 
the orbital period, log g and chemical abundances of the stars 
that pass these selection criteria with the observations. For each 
of the 15 observed stars in our sample the best-fitting model is 
found with the following iterative procedure. 


mass dredged up during the AGB phase, M DU , tot , as a function of 
the initial mass of the star (top and bottom panels, respectively). 

Because the protons in the PMZ are processed to produce 
13 C, the parameter Mpmz plays the same role in the chemical 
evolution of the AGB star as the efficency of the 13 C pocket dis- 


1. We first select the modelled binary systems with the right 
orbital period. Because the observed errors on the periods 
are generally small (or in some cases this information is not 
provided in the literature) the working error that we adopt 
coincides with the spatial resolution in our grid of models, 
which corresponds approximately to an uncertainty in the 


cussed by e.g. Straniero et al. 

(1995 

), Busso et al. ( 

1999) and, 2. Among the modelled stars that belong to binary systems 

more recently, Bisterzo et al. 

(2010 

). The main effects of in- within the correct range of orbital periods, we select those 


creasing Mpmz in our model, other conditions being equal, are 
to bring more 5-elements to the stellar surface and to produce a 
distribution of 5-elements increasingly weighted towards lead. In 
Fig. [3]we show the element abundances in a 2 M 0 star at the end 
of the AGB computed with binary c for four different masses 
of the PMZ. While the abundances of light elements differ by 
about 0.5 dex at most (Na, Mg, P) and generally remain con¬ 
stant, the abundances of 5-elements vary significantly between 
the minimum to the maximum size of the PMZ, M pm z = 0 M 0 
and Mpmz = 0.004 M 0 , respectively. In Fig. [3] we show the re¬ 
sults of the Stromlo model for comparison (plus signs): the abun¬ 
dances predicted by the two codes agree to within 0.1 dex. We 
note that below about 2M 0 the surface abundances are much less 
sensitive to Mpmz and generally only differ by few 0.1 dex, be¬ 
cause for low masses proton-ingestion episodes occur in the de¬ 
tailed models that cause 5-elements production independently of 


Mpmz (we refer to Lugaro et al. 2012 for details). 


4.3. Method 

As previously mentioned, the aim of our work is to find for each 
star in our observed sample a model that reproduces the mea¬ 
sured abundances. To this end, we generate a grid of evolution 
models of A binary stars distributed in the M\ - M 2 - log 10 a - 
Mpmz parameter space, where M\ 2 are the initial masses of the 
primary and of the secondary star, respectively, a is the initial or¬ 
bital separation of the system and M pmz is the mass of the partial 
mixing zone of any star that undergoes AGB evolution. The grid 
resolution is equal to A - Ami x Am 2 x N a x Apmz, where we 


stars that reproduce the observed surface gravity within one 
cr g , where cr g is the observed error on log g. The value of 
log g gives an indication of the phase of evolution of a star. 
The constraint on log g guarantees that we reproduce the evo¬ 
lutionary status of the observed stars. 

3. For each of the modelled stars that have passed the selection 
criteria on P or b and log g we calculate the x 2 for the model 
fit to the abundances with the following equation: 


X 


-Z 


(A^obs A/.mod)' 


07 


(4) 


where A, i0 b S and A/ mo d are the abundances of element i (de¬ 
fined as in Eq. |T| in the observed and modelled star, respec¬ 
tively, and cr, is the observed uncertainty. We are mostly in¬ 
terested in studying the nucleosynthesis that occurs during 
the AGB phase, so for this reason we compute x 2 taking into 
account only the elements that are produced or destroyed 
by AGB stars. This includes the light elements C, N, O, F, 
Ne, Na, Mg and all the heavy neutron-capture elements with 
atomic number in the range 31 to 82. In this range are in¬ 
cluded the peak of light-5 elements around atomic number 
40 (strontium, yttrium and zirconium), the peak of heavy- 
5 elements around atomic number 56 (barium, lanthanum 
and cerium) and the lead peak. In AGB stars of mass below 
~ 3 M 0 nitrogen is produced by the CN cycle when protons 
from the envelope are mixed into regions of the star where 
carbon is abundant. The exact amount of nitrogen depends 
on the extent of this mixing process that is uncertain and 
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Fig. 3. Element abundances relative to iron predicted for a star of mass M, = 2M 0 with Mpmz = 0, 2 x 10 4 , 2 x 10 3 , 4 x 10 3 M 0 
(respectively: dotted, dot-dashed, dashed, solid line). Plus signs represent the value predicted by the Stromlo model for the same 
MpMZ- 


is not included in our AGB nucleosynthesis model; how¬ 
ever, the total amount of carbon and nitrogen is conserved 
and therefore when both the elements are observed we fit 
our models to the combined abundance C+N. In our study 
we do not consider the elements with atomic number be¬ 
tween 13 (aluminium) and 30 (zinc). The abundance of alu¬ 
minium is not affected by low-mass AGB stars ( M, < 3M 0 ), 
which are the dominant site of the s-process, and its pro¬ 
duction mostly occurs in intermediate-mass AGB stars with 
hot-bottom burning, M, > 3 M Q (| Ventura & D’Antona 2009 


Karakas et al. 2012[>. Elements heavier than aluminium up 


to zinc are not in general produced by AGB stars of any 
mass range (Karakas et al. 2009 1 Cristallo et al. 2011 1 . The 
abundances of these elements in very metal-poor stars are 
expected to be consistent with the chemical composition of 
the gas cloud from which the stars were formed and to be 

reproduced by our adopted set of initial abundances, that are 

based on the results of galactic chemical evolution models 
at metallicity Z ~ 10 * * * 4 . The discrepancies between the ob¬ 
served abundances of elements between aluminium and zinc 
and our set of initial abundances are discussed more in detail 
in Paper II. 

4. The model that gives the best match to the observed abun¬ 
dances is determined by the minimum value of x 2 , X^in 
estimate of how well this model reproduces the observed 
abundances is given by the reduced^ 2 , i.e. xl = A 2 /A where 
v is the number of degrees of freedom. This is calculated as 


v — (V 0 bs (n [it n c ) — (V 0 bs 3 , 


(5) 


where (V 0 b s is the number of observables, i.e. the number of 
elements that are used to calculate ^ in ; n b t = 5 is the num¬ 
ber of the fitted parameters, i.e. M\, Mi, M pmz , log g, P or b; 
and n c - 2 is the number of observational constraints that 
are not directly involved in the calculation of x 2 although 
they limit the range of acceptable models, i.e. logg and the 
observed orbital period. Upon visual inspection models with 
/v < 3 appear to fit the observed abundances well, while 


above this threshold the observations are not reproduced 
well. The reduced^ 2 should not be used for a goodness-of-fit 
statistical test, because the uncertainties associated with the 
observed abundances generally only convey the uncertainties 
of the method adopted to determine the abundances, while 
systematic errors are not taken into account, which may be 
caused for example by the errors in the estimates of the tem¬ 
peratures and gravities of the stars (as we mentioned, for ex¬ 
ample, for star CS22948-027). 

5. We determine how well-constrained our best model is by 
calculating the confidence intervals of its input parameters. 
Mix Mix Pi and Mpmz- On our grid of binary models and 
associated x 2 values, we fix one of the input parameters p 
and we find the minimum x 2 with respect to the other input 
parameters. This procedure defines a function that associates 
a^ 2 to each grid value of the parameter p. The minimum of 
this function is equal to the x 2 of the best fit, X^in- ^ confi¬ 
dence region is defined as an interval of p within which the 
difference^ 2 - Amin' s below a certain threshold. If the mea¬ 
surement errors are Gaussian and our model reproduces the 
data, then the probability distribution of a 2 is also Gaussian. 
In this case, the threshold A^ 2 = 1 corresponds to the con¬ 
fidence interval of 68.3% probability that the actual p is in 
this interval. Similarly, the thresholds Ax 2 = 4 and A^ 2 = 9 
correspond to the confidence intervals of 95.4% and 99.7% 
probability, respectively. However, because systematic errors 
are generally not taken into account in the estimated obser¬ 
vational errors, the thresholds A;^ 2 = 1,4,9 should not be 
used to calculate the theoretical Gaussian probabilities. In 
the analysis of our sample stars we note that models with 
X 2 below the threshold Ax 2 = 4 have surface abundances 
that are hard to distinguish by eye from the best-fit values, 
whereas higher x~ are found for models clearly distinct and 
worse than the best model. Therefore, to determine the con¬ 
fidence intervals of our input parameters we adopt A^ 2 = 4. 
In case for one parameter no model other than the best fit 
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satisfies the condition Ay 2 < 4, the confidence interval is 
assumed to be half of the grid resolution for that parameter. 


In points 1 - 4 of our procedure we do not include any con¬ 
straint on r eff , which has a strong dependence on the metallicity. 
The observed metallicities of the stars in our sample vary by 
up to a factor of five. On the other hand, in our model we keep 
the metallicity constant because the nucleosynthesis is valid for 
Z = 10 4 and hence we do not expect to be able to reproduce T e ff. 

In our study we do not take into account that most neutron- 
capture elements are not purely produced by the s-process 
but also have an r-process component (Arlandini et al.|[l999| 


Bisterzo et al. 2011 1 . The origin of the r-process elements in 


meta 

-poor stars is still unclear (Sneden et al. 1994; lonsell 

et al. 

2006 

Sneden et al. 2008; [Lugaro et al. 2009; [Masseron 

et al. 

2010 

and a thorough analysis of this issue is beyond the 


scope of this paper. 


5. Results 

We determine the best-fitting models to the observed abundances 
of the 15 stars in our sample. In Tables [2] and [3] we summarise 
the results obtained with model sets A and B, respectively, as 
follows. 


• Columns 2-6: the initial masses of the modelled stars, 
the mass of the PMZ, the initial orbital period and, for bi¬ 
nary systems that undergo common-envelope evolution, the 
adopted value of the common-envelope efficiency. 

• Columns 7 - 12: the results of the fit, namely the period 
of the modelled binary, Pf, at the moment when the mod¬ 
elled secondary star best reproduces log g and the observed 
abundances, the surface gravity, the effective temperature, 
the amount of mass accreted by the secondary AM acc , the 
minimum value of x 2 , X 2 mn > the number of degrees of free¬ 
dom of the fit, v, and the reduced A 2 ’Am in / V - 


In the next section we describe the main results of our study, 
while in Sect. l5.2l and l5.3l we discuss in detail the models of the 
r-normal CEMP-s star CS22942-019 and the CEMP-s/r star 
CS29497-030 as an example of our analysis. In Sect. 5.4 we 


discuss the method adopted to model the CEMP-s binary star 
CS22964-161A,B. The plots of the best-fitting models of the re¬ 
maining 11 stars of our sample are shown in Appendix [A| 


5. 1. Main results 

The observed orbital periods provide strong constraints to the 
initial orbital periods in our models and consequently to the ini¬ 
tial primary masses. In our model set A (with a spherically sym¬ 
metric wind), the systems typically widen in response to mass 
loss. Hence, the final orbital periods reproduce the observations 
only if the initial periods are relatively short (less than a few 
thousand days). Consequently, the initial primary masses are rel¬ 
atively low (generally M \j < 1.1 M 0 ), because the primary stars 
have to evolve up to the AGB without overfilling their Roche 
lobes, otherwise the binary systems start a common-envelope 
phase without mass accretion. On the other hand, with model 
set B (efficient angular momentum loss) the systems shrink in 
response to mass loss, therefore longer initial orbital periods are 
possible and consequently higher primary masses. 

Six stars in our sample have periods shorter than 500 days. 
For these systems our model set B predicts that initially the pe¬ 
riods are up to a hundred times longer. During the AGB phase 


the primary star loses about 60 - 70% of its envelope mass in 
the wind before filling its Roche lobe. The envelope mass that is 
left forms a common envelope, the ejection of which requires the 
system to lose angular momentum and shrink to the observed or¬ 
bital period. In star HD224959 (Fig. |A.8| > the common-envelope 
phase is preceded by a short phase of stable Roche-lobe over¬ 
flow, during which 0.1 M 0 are accreted by the secondary star. 
Star HE0024—2523 (Fig. | A.9 [ > has an observed orbital period of 
3.14 days. In this close orbit all primary stars of our grid overfill 
the Roche lobe while ascending the red giant branch. An ini¬ 
tially wider orbit is necessary for the primary star to evolve to 
the AGB phase. During this phase some material is transferred 
onto the secondary star before the system enters in a common 
envelope that is ejected with very low efficiency (a cm = 0.07 
or 0.03 for model sets A and B, respectively). An alternative in¬ 
terpretation that does not require such a low common-envelope 


efficiency is discussed in Sect. 6.1 


Eight stars in our sample have passed the main-sequence 
turnoff and have low surface gravity. As a consequence, large 
amounts of accreted material are required to reduce the dilution 
caused by the first dredge-up and reproduce the enriched abun¬ 
dances of carbon and ,v-elements. At periods shorter than about 
5,000 days small amounts of mass are accreted with our model 
set A, while in model set B (with enhanced BHL accretion effi¬ 
ciency) the accreted mass is up to a few tenths of a solar mass. As 
a consequence of the constraints on the period, M \and AM acc , 
the best-fitting models to the abundances of ten stars are found 
with model set B. Model set A predicts a significantly better fit 
than set B only for star BD+04°2466, which is the longest-period 
binary in our sample and exhibits relatively small enhancements 
([El/Fe] < 2 for all observed elements. Fig. |A.1| >. Hence, a rel¬ 
atively small accreted mass (AM acc < O.1M 0 ) is sufficient to 
reproduce the observed abundances. 

We obtain X 2 m [J v - 3 only for six stars. The confidence in¬ 
tervals of the input parameters of these systems are summarised 
in Table |d] Multiple local minima occur in the x 1 distributions 
of the primary masses and initial separation of stars BD04°2466 
and HD201626. These minima are found because different com¬ 
binations of initial parameters result in model stars with similar 
surface abundances and hence similar x 2 - A few examples of 
the confidence intervals determined in our study are discussed in 
Sect. 5.2 and in Appendix [5] 


5.2. Example 1: the r-normal CEMP-s star CS22942-019 

In Fig. [4] we show the observed and modelled abundances of 
star CS22942-019. Black points with solid error bars show the 
abundances of the elements taken into account in Eq. Q to de¬ 
termine the best fit. The other observed elements are shown in 
grey. Our best model (red solid line) is found with model set B, 
in which a primary star of mass M\ \ = 1.8 M 0 transfers material 
to a 0.54 M q companion (Table [3]). In the lower panel of Fig. [4] 
we show the residuals calculated as the difference between the 
observed and the modelled abundance of each element. Our best 
model is found combining three ingredients: an initially wide 
separation, so that the primary star does not fill its Roche lobe, 
high mass accretion rate, and efficient angular momentum loss, 
which is required to shrink the initially wide system to the ob¬ 
served period. 

In our best model computed with our default model set A 
(green dotted line in Fig. [4} the binary system widens instead 
of shrinking as in model set B. Consequently, the initial or¬ 
bital separation needs to be shorter (P; ~ 2100 days) to avoid 
filling the Roche lobe and therefore a lower M\\ is necessary 
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Table 2. Physical parameters of the modelled stars computed with model set A with WRLOF wind-accretion rate and spherically 
symmetric wind. 


ID 

M u 

M 2 ,; 

Mpmz 

P ; 

ffCE 

Pi 

l°gio g 

7eff 

AA/acc 

y 2 . 

'A. min 

V 

xiJy 

BD+04°2466 

1.1 

0.76 

2 x 10 3 

3.19 x 10 3 


4.28 x 10 3 

1.61 

4600 

0.09 

13.4 

8 

1.7 

CS22942-019 

1.4 

0.76 

4 x 10' 3 

2.10 x 10 3 


2.87 x 10 3 

2.60 

5000 

0.10 

74.3 

10 

7.4 

CS22948-027 

0.9 

0.81 

2 x 10' 3 

4.19 x 10 2 


4.76 x 10 2 

2.20 

4900 

0.12 

82.8 

10 

8.3 

CS22956-028 

0.9 

0.84 

2 x 10" 3 

1.17 x 10 3 


1.37 x 10 3 

3.82 

6900 

0.10* 

50.9 

2 

25.4 

CS22964-161A 

1.6 

0.79 

2 x 10" 3 

2.74 x 10 5 


5.60 x 10 5 

3.50 

5800 

0.04 

19.5 

8 

2.5 

CS22964-161B 

1.6 

0.71 

2 x 10" 3 

2.74 x 10 5 


5.60 x 10 5 

4.48 

6200 

0.02 

11.8 

9 

1.3 

CS29497-030 

0.9 

0.84 

2 x 10' 3 

3.37 x 10 2 


4.04 x 10 2 

3.63 

6500 

0.14* 

145.6 

17 

8.6 

CS29497-034 

1.9 

0.76 

2 x 10' 3 

2.67 x 10 3 


3.35 x 10 3 

1.80 

4700 

0.11 

117.6 

10 

11.8 

CS29509-027 

2.9 

0.71 

2 x 10" 3 

3.09 x 10 3 

1.0 

1.65 x 10 2 

4.14 

7000 

0.08 

4.6 

1 

4.6 

HD198269 

0.9 

0.79 

3 x 10~ 3 

1.19 x 10 3 


1.43 x 10 3 

1.55 

4600 

0.06 

32.4 

6 

5.4 

HD201626 

0.9 

0.81 

2 x 10' 3 

4.19 x 10 2 


4.75 x 10 2 

2.23 

4900 

0.12 

30.1 

6 

5.0 

HD224959 

0.9 

0.79 

3 x 10" 3 

1.19 x 10 3 


1.42 x 10 3 

2.05 

4800 

0.06 

340.1 

9 

37.8 

HE0024-2523 

1.9 

0.71 

2 x 10' 3 

2.57 x 10 3 

0.07 

3.33 

4.40 

6500 

0.06 

203.6 

5 

40.7 

HE0507-1430 

0.9 

0.81 

0 

4.19 x 10 2 


5.16 x 10 2 

0.70 

4300 

0.12 

26.5 

0 

- 

LP625-44 

1.9 

0.76 

4 x 10" 3 

2.67 x 10 3 


4.03 x 10 3 

2.95 

5100 

0.11 

374.0 

18 

20.8 


Table 3. Physical parameters of the modelled stars computed with model set B in which we adopt an efficient Bondi-Hoyle-Lyttleton 
wind-accretion rate (ctbhl = 10) and efficient angular-momentum loss. 


ID 

M u 

M 2 

Mpmz 

Pi 

CCE 

Pi 

lo gio£ 

Tee 

AM acc 

y 2 . 

'A. min 

V 

xLJy 

BD+04°2466 

0.9 

0.79 

3 x 10~ 3 

4.74 x 10 3 


4.25 x TO 3 

1.70 

4700 

0.05 

25.7 

8 

3.2 

CS22942-019 

1.8 

0.54 

2 x 10 4 

1.67 x 10 4 


3.23 x 10 3 

2.26 

4900 

0.29 

21.0 

10 

2.1 

CS22948-027 

1.5 

0.61 

1 x 10' 3 

6.72 x 10 3 

1.0 

3.69 x 10 2 

2.20 

4900 

0.27 

31.9 

10 

3.2 

CS22956-028 

0.9 

0.84 

2 x 10' 3 

1.17 x 10 3 


1.05 x 10 3 

3.90 

7100 

0.10* 

56.3 

2 

28.2 

CS22964-161A 

1.6 

0.79 

2 x 10' 3 

1.30 x 10 5 


1.56 x 10 5 

3.50 

5800 

0.04 

19.6 

8 

2.5 

CS22964-161B 

1.6 

0.71 

2 x 10' 3 

1.30 x 10 5 


1.56 x 10 5 

4.48 

6200 

0.02 

11.8 

10 

1.3 

CS29497-030 

1.5 

0.56 

2 x 10" 3 

7.17 x 10 3 

1.0 

3.00 x 10 2 

4.34 

6700 

0.22* 

116.0 

17 

6.8 

CS29497-034 

1.6 

0.56 

2 x 10' 3 

1.18 x 10 4 


3.35 x 10 3 

1.77 

4700 

0.31 

34.1 

10 

3.4 

CS29509-027 

2.9 

0.59 

5 x 10~ 4 

2.62 x 10 4 

1.0 

1.67 x 10 2 

4.28 

6900 

0.19 

5.0 

1 

5.0 

HD198269 

1.0 

0.81 

3 x 10" 3 

1.62 x 10 3 


1.17 x 10 3 

1.05 

4400 

0.13 

18.0 

6 

3.0 

HD201626 

2.6 

0.64 

4 x 10" 3 

2.72 x 10 4 

1.0 

3.55 x 10 2 

2.50 

5000 

0.22 

13.6 

6 

2.3 

HD224959 

1.2 

0.61 

4 x 10' 3 

3.24 x 10 3 


1.41 x 10 3 

2.12 

4800 

0.34 

72.2 

9 

8.0 

HE0024-2523 

1.1 

0.69 

0 

2.15 x 10 3 

0.03 

3.48 

4.40 

6600 

0.11 

50.0 

5 

10.0 

HE0507-1430 

1.8 

0.61 

0 

1.12 x 10 4 

1.0 

3.58 x 10 2 

0.70 

4300 

0.28 

7.6 

0 

- 

LP625-44 

1.8 

0.56 

4 x 10' 3 

1.66 x 10 4 


3.65 x 10 3 

2.94 

5100 

0.30 

194.4 

18 

10.8 


Notes. All the masses are expressed in units of M 0 , periods are in days, r eff in K, and g is in units of cm s~ 2 ; v is the number of degree of freedom 
of the fit; a CE is shown for modelled systems that experience a common-envelope phase. 

(*) inefficient thermohaline mixing, the accreted material stays on the stellar surface. 


Table 4. Confidence intervals of the input parameters of our model stars with x 2 ! v - 3. 


ID 

model set 

xiJy 

best 

M u /M 0 

min 

max 

best 

M 2 ,i/M 0 

min 

max 

Mpmz/( 10 
best min 

J M 0 ) 

max 

Pi/(10 3 days) 
best min max 

BD+04°2466 

A 

1.7 

1.10 

0.95 

1.15 

0.76 

0.73 

0.85 

2.0 

0.05 

4.0 

3.19 

2.32 

3.92 





1.45 

1.65 










CS22942-019 

B 

2.1 

1.80 

1.75 

1.85 

0.54 

0.53 

0.58 

0.2 

0.0 

0.4 

16.7 

14.1 

19.8 

CS22964-161A 

B 

2.5 

1.60 

0.95 

1.65 

0.79 

0.75 

0.85 

2.0 

0.58 

4.0 

130 

103 

232 

CS22964-161B 

B 

1.3 

1.60 

0.95 

1.65 

0.71 

0.68 

0.85 

2.0 

0.58 

4.0 

130 

103 

232 

HD 198269 

B 

3.0 

1.00 

0.95 

1.05 

0.81 

0.73 

0.83 

3.0 

0.80 

4.0 

1.62 

1.36 

1.93 

HD201626 

B 

2.3 

2.60 

1.95 

2.95 

0.64 

0.58 

0.65 

4.0 

0.80 

4.0 

27.2 

15.2 

36.3 





1.55 

1.65 








6.98 

9.86 


(Mi,; = 1.2 - 1.4 M 0 ). With these initial parameters our WRLOF 
prescription is not very efficient: the secondary star accretes a 
smaller amount of material compared to our best model, approx¬ 
imately A/Mjcc =0.1 M 0 , and therefore it needs to be initially 
more massive, M 2j - 0.76 M 0 , to reproduce the observed log g. 
The accreted material is more diluted in the secondary star and 
we underestimate the observed abundance of carbon, the light-s 
elements and barium even when assuming the largest PMZ avail¬ 


able in our model, Mpmz = 4x 10 Consequently the fit is much 
poorer, X^ m = 74. 

The choice of a relatively massive primary star in model 
set B (Mi,; = 1.8 M 0 ) implies that the abundances of almost 
all the elements are reproduced with Mpmz = 2 X 10 4 . A pri¬ 
mary star more massive than 1.8 M 0 overproduces sodium (even 
with Mpmz = 0) and typically produces excessively abundant 
heavy-s elements, in contrast with the observations, which give 
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Fig. 4. Best fitting model to star CS22942-019. Points with error bars: observed abundances, in black the elements used to determine 
the best fit and in grey the other elements. Red solid line: the best-fitting model, found with model set B and My; = 1.8 M 0 , 
M 2 ,i = 0.54 M 0 and Mpmz = 2 x 10 4 M 0 . Blue dashed line: alternative fit with model set B and My; - 1.3 M 0 , M 2j ; = 0.61 M 0 and 
Mpmz = 2 X 10 3 M 0 (see text). Green dotted line: best fit adopting model set A. Lower panel: the residuals of the three models are 
shown as red plus signs with error bars, blue triangles, and green dots, respectively. 





Fig. 5. One-dimensional confidence intervals of the input param¬ 
eters for model star CS22942-019. Panels a-d show the initial 
primary and secondary masses. My; and M 2> i, the PMZ mass, 
Mpmz- and the initial period, / J ,, respectively. Long-dashed lines 
indicate the thresholds A^ 1 = 1, 4, 9. 


[hs/ls] » 0. Lugaro et al. (2012) consider M pmz = 2x 10~ 3 M 0 
as the standard PMZ mass and most models in Tables [2] and Rl 
have a relatively large PMZ (typically few 10 3 M 0 ). In Fig. [4] 


we compare our best fit with the best model that adopts Mpmz = 
2 x 10 3 M 0 (blue dashed line). The initial parameters of this 
model are My; = 1.3 M 0 , M 2 ,i = 0.61 M 0 , P ; = 6.28 x 10 3 days, 
and the other assumptions are the same as in the best model. 
With this model the light-5 1 elements are underestimated while 
most of the heavy-.v elements are overestimated and as a result 
X 1 = 54. 

In Fig. [5] we show the confidence intervals of the input 
parameters of model star CS22942-019 computed with set B. 
Panels a-d show the one-dimensional confidence intervals of the 
initial primary and secondary masses, PMZ mass and initial or¬ 
bital period, respectively. No model is found with Ax 1 < 9 if 
Myi or P[ differ from the best-fitting values and few models ful¬ 
fil the condition A_y 2 < 4 for a variation of M 2j or Mpmz- These 
results do not imply that the initial parameters of the progeni¬ 
tor system of CS22942-019 are determined with small uncer¬ 
tainty. The narrow confidence intervals that we derive (Table [4} 
are a consequence of the constraint put by the observed period. 
With a 1.8M 0 primary star we find only one initial period in our 
grid of models for which a binary system evolves into a system 
with the observed orbital period. If a different initial primary 
mass is selected, no combination of the other three parameters 
is found that reproduces the measured period and the observed 
abundances with A^ 1 < 4. In Paper II we show that the confi¬ 
dence intervals of Mpmz. Pi and M 2 ,; are larger if we release the 
period constraint. 

5.3. Example 2: the CEMP-s/r star CS29497-030 

The observed period of star CS29497-030 in Fig. [6] is rather 
short, P or b = 342 days. In our model set B the binary system 
is initially wider (about 7200 days). The primary star (My; = 
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Fig. 6. Points with error bars: as Fig. |d] for star CS29497-034. Red solid line: best-fitting model found with model set B and 
inefficient mixing of the material accreted by the secondary star. Blue dashed line: same input parameters as the best model but 
assuming efficient thermohaline mixing. In the right panel the combined abundance of carbon and nitrogen is shown. 


1.5M G ) ascends the AGB and loses approximately 0.53M 0 in the 
wind, which is partly accreted by the secondary star (AM acc = 
0.22M 0 ). Subsequently, when its residual envelope mass is about 
0.25M q , the primary star fills its Roche lobe and the system en¬ 
ters a common-envelope phase. The ejection of the common en¬ 
velope, modelled with acE = 1, shrinks the system to the period 
Pf = 300 days that reproduces the observed value within the 
uncertainty of our grid. The common-envelope phase is required 
because the wind of the primary star does not carry away enough 
angular momentum even with our model set B. 

With our model set A we find a binary system with almost 
equal initial primary and secondary masses (Mf; = 0.9M 0 , 
Mix = 0.84M o ) and short initial period (P, = 337 days) that ex¬ 
pands as an effect of mass loss and does not undergo a common- 
envelope phase. The amount of mass transferred is small com¬ 
pared to model set B (A M acc = 0.14M 0 ) and the fit is signifi¬ 
cantly worse 0f^ in ~ 146, whereas x^ mn = 116 for model set B). 

Star CS29497-030 is strongly enhanced in carbon, heavy- 
s elements, /--process elements and lead. On the other hand, 
sodium, magnesium and the light- s elements are only weakly 
or mildly enhanced (between 0.4 and 1.5 dex). Similar abun¬ 
dance distributions are also observed in the other CEMP-s/r 
stars of our sample (CS22948-027, CS29497-034, HD224959 
and LP625-44). Our best model (red solid line in Fig. [6} re¬ 
produces the abundances of barium, most /--elements and lead 
within the observational uncertainty but it overestimates the 
abundance of sodium (by approximately 1.5 dex), light-5 1 ele¬ 
ments, lanthanum and cerium (by 0.3 - 1 dex). Our best model 
is computed assuming that non-convective mixing mechanisms 
of the accreted material, such as thermohaline mixing, are ineffi¬ 
cient. Consequently, the transferred material remains on the sur¬ 
face of the secondary star because the relatively high observed 
gravity, log I0 (g/cms~ 2 ) = 4.0, indicates that CS29497-030 is 


still on the main sequence and therefore it has not yet undergone 
the first dredge-up. 

The blue dashed line in Fig. [6] shows the abundances pre¬ 
dicted by a model with the same input parameters as in the 
best model and efficient thermohaline mixing. This model re¬ 
produces the abundances of C+N, Mg, light-s elements. La, Ce 
and Nd better than the best-fitting model, but the x 2 is higher 
(X 2 = 167) because it fails to reproduce the eight elements be¬ 
tween europium and lead. Except for lead, which is off by less 
then 2cr, these elements are mostly produced by the /'-process. 
If we exclude the /'-elements from the calculation of the fit we 
find the same initial parameters as in the best model. The re¬ 
duced^ 2 is still high Of 2 /v = 7.4) because it is dominated by 
the discrepancy in the abundance of sodium, which is overesti¬ 
mated by approximately 1 dex. Because both models in Fig.[6]do 
not reproduce the observed abundances well Ofmi J v > 3) we do 
not have strong constraints on the efficiency of non-convective 
mixing processes in CS29497-030. 


5.4. Example 3: the CEMP-s binary star CS22964-161A,ES 


The binary system CS22964-161A,B represents a special case, 
because both stars are enriched in carbon and ,y-elemcnts and 
have relatively high surface gravities, consistent with stars that 
have recently passed the main-sequence turnoff (star A, Fig.[7j0 
or are still on the main sequence (star B, Fig. [ 7 J 5 ). The simplest 
interpretation is that the binary system seen today was the recip¬ 
ient of the material lost in the past from a third star in its AGB 
phase, in a hierarchical triple system (Thompson et al. 2008). 
In this hypothesis the observed orbital period (P or t, = 252 days) 
does not correspond to the final period of our models, Pf, which 
instead is the period of the conjectured, unseen third star. 

Our binary_c code is suited to study the evolution of bi¬ 
nary systems and not triple systems, hence to determine the best- 
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Fig. 7. Best-fitting models and observed abundances of stars CS22964-161A (left) and CS22964—161B (right). Symbols are the 
same as in Fig. [4] 


fitting model to the observed abundances we use the following 
method. We find the models that best reproduce the abundances 
of stars A and B separately, with no restriction on the orbital 
period but with the constraint that the primary mass has to be 
the same for the two stars. We determine the initial masses Mu, 
Mrmz, M a .j and Mb. j. We find that stars A and B accrete 0.04M o 
and 0.02M e , respectively, and the total x^ in is X^ mn A+B — 31.4 


(with model set B). Subsequently we study the evolution of a bi¬ 
nary system in which the initial primary and secondary masses 
are Mg and M^x = My, + Mb.;, respectively. In this way we 
mimic the effect of a triple system where the primary star in a 
wide orbit “sees” the close inner binary as a single object with 
mass equal to the sum of the two components. We impose that 
the primary star transfers to the companion an amount of mass 
equal to M acc = 0.06M G and we determine the initial orbital 
period. In our best model star A is initially 10% more massive 
than star B (Max = O.79M 0 and Mb.; = 0.71M o ). The final 
mass ratio q = Mb/M a is approximately 0.90, similar to the 
value estimated from the observations (q = 0.88, |Thompson 
et al., 2008)>. A choice of two stars with initially equal masses 


Ma = Mg = 0.76 M Q that accrete the same amount of material 
AM acc ~ 0.04 M 0 leads to a model that is similar to our best 
model (x\ +B =33). 

Our best model reproduces the abundances of all the ele¬ 
ments in both stars, except for sodium, which is overproduced 
in star B, and yttrium, which is slightly overestimated in both 
stars. We find equally good fits for both model sets, with the 
same parameters except for the initial and final periods of the 
third star. 


6. Discussion 

6.1. Constraints on mass transfer and binary evolution 

Our model set B (with enhanced Bondi-Hoyle-Lyttleton wind- 
accretion efficiency and efficient orbital angular momentum loss) 
provides the best fit to the observed abundances in ten stars and 
predicts results very similar to the best fit in four of the remain¬ 
ing five systems. This model set is not necessarily realistic be¬ 
cause we assume an arbitrarily high efficiency of wind mass ac¬ 
cretion. However, these results indicate that the accretion to the 
secondary star of large amounts of material is necessary in at 
least ten stars of our sample, otherwise we find a poor fit to the 
observations. Eight of these stars have low surface gravity, which 


implies that they evolved off the main sequence and have under¬ 
gone the first dredge-up. Consequently, the accreted material is 
diluted throughout the envelope of the recipient star, regardless 
of our assumptions about the efficiency of thermohaline mix¬ 
ing. Large amounts of transferred mass are therefore required 
to reduce the dilution and reproduce the large enhancements of 
carbon and .v-elements. 


Our model set A (with WRLOF wind-accretion efficiency 
and spherically symmetric wind) is disfavoured by the require¬ 
ment of large amounts of mass accretion combined with the con¬ 
straint on the orbital period. Because with a spherically symmet¬ 
ric wind the orbit typically expands in response to mass loss, 
the modelled binary stars need to be initially close. For orbital 
periods shorter than a few thousand days even a low-mass pri¬ 
mary star transfers a small amount of material to the compan¬ 
ion in our WRLOF model. Consequently, with model set A only 
stars BD+04°2466 and CS22964-161A,B are well reproduced 
because the orbital periods are long and the observed carbon and 
s-elements are not strongly enhanced. In the future more realistic 
models of wind mass transfer are required to take into account 
more accurately the accretion process in close binary stars. 


In six systems a common-envelope phase is necessary to 
shrink the orbit to the observed period (assuming model set 
B). The best models of five of these systems adopt the default 
common-envelope efficiency qtce = 1- Star HE0024—2523 (Fig. 


A.91 has an observed period of three days that can be reproduced 


only by assuming a very inefficient ejection process, acE = 0.03. 
An alternative interpretation is that HE0024-2523 was initially 
part of a hierarchical triple system in which an intermediate- 
mass primary star was in a wide orbit around two low-mass 
stars in a close binary. In this scenario the inner binary enters 
in a common envelope after being polluted by the primary star. 
Therefore, to reproduce the observed period an inefficient pro¬ 
cess of common-envelope ejection is not necessary. In an alter¬ 
native triple scenario the unseen companion of HE0024—2523 
is a low-mass main-sequence star and the inner binary was ini¬ 
tially formed with the observed period of three days, in which 
case a common-envelope phase is not required. In summary, the 
acE parameter is not well constrained and we do not find strong 
evidence to reject the default value oqe = L 
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6.2. Comparison between modelled and observed 
abundances 

In some systems our model predictions do not reproduce the ob¬ 
served abundances even when we adopt a fine-tuned model of 
wind mass transfer. The abundance of nitrogen, for example, is 
determined in eight stars of our sample and it is well reproduced 
by our model only in two systems (HD201626 and HE0507- 
1430 in Figs. A.7 and A. 10 respectively). However, as we noted 
in Sect. 4.3 the exact amount of nitrogen produced by the CN cy¬ 
cle in AGB stars is uncertain. If some extra amount of carbon is 
converted to nitrogen in the model we reduce the discrepancies 
and in six out of eight stars the abundance of C+N is reproduced 
within the observational uncertainty. 

The abundance of oxygen does not vary by more than ap¬ 
proximately 1 dex in the models of low-mass AGB stars and 
the models typically underestimate the observations. This dis¬ 
crepancy may indicate that some oxygen is mixed to the surface 
by the third dredge-up or that the initial abundance of oxygen 
adopted in our models is low compared to the average value ob¬ 
served in our sample of very metal-poor stars, as we will discuss 
in Paper II. 

In many stars the observed abundances of sodium are low 
compared to the predictions (e.g. CS29497-030 and CS22964- 
161B). Sodium is produced in the intershell region by proton 
capture on neon seeds, therefore the sodium abundance grows 
very rapidly with increasing PMZ mass (see also the discus- 

and references 


Lugaro et al. 


2012 


sion about 23 Na in K10 and 
therein). A better match would generally require a lower mass 
of the PMZ, but this is hard to reconcile with the large en¬ 
hancements of ,v-elements that need a massive PMZ to be repro¬ 
duced. The abundances of sodium in low-metallicity stars show a 
large dispersion which is difficult to reproduce in galactic chem¬ 
ical evolution models (Kobayashi et al.[ |2011| l. However, our 
AGB models of mass above « 1.2 M Q produce large amounts 
of sodium and a factor of ten difference in its initial abundance 
accounts for no more than 0.1 dex in the final [Na/Fe]. Therefore 
we consider it more likely that the discrepancy in sodium is re¬ 
lated to the large uncertainties in the numerical treatment of the 


PMZ in the detailed models (e.g., [Goriely & Mowlavi 
Lugaro et al., 2004j >. 


2000 


In CS22964-161A,B the observed abundances are generally 
well reproduced, but yttrium is several tenths of a dex lower than 
strontium and zirconium. Such a distribution of strontium, yt¬ 
trium and zirconium is observed in other CEMP-5 stars, as we 
will discuss in Paper II, but it is not replicable by our mod¬ 
els because these elements are always produced in comparable 
amounts. The observed abundance of zirconium is low compared 
to the models of the stars HD198269, HD201626 and HD224959 
(Figs. A.6f A.8 >. This discrepancy may be related with a sys¬ 
tematic effect in the observations by |Van Eck et al.| ( |2003| ), be¬ 
cause previous measurements from spectra at lower resolution 
by Vanture (1992a|b I are consistent with the models. 


6.3. Abundances in CEMP-s/r stars 


The CEMP-.v/r stars in our sample, CS22948-027 (Fig. A.2i, 
CS29497-030 (Fig. [6), CS2949 7-034 (Fig. \AA\ . HD224959 
(Fig. A.8 i and LP625-44 (Fig. A. 1 1 [>, show highly enriched 


abundances of elements that are mostly produced in the r- 
process, such as europium. These elements are typically pro¬ 
duced in small amounts in our AGB nucleosynthesis model. 
Consequently, to reproduce the large enhancements of the r- 
elements our models often overestimate the abundances of other 


elements, such as carbon, sodium and light-5 elements. This is¬ 
sue illustrates the general problem of reproducing element-to- 
element ratios. More specifically, in some systems a high value 
of Mpmz is necessary to match one element whereas a lower 
Mpmz is sufficient to reproduce another element. In LP625^44, 
for example, all the elements in the heavy-5 peak are underes¬ 
timated and would need a higher Mpmz (or a higher M\) to be 
reproduced. On the contrary light elements, light-5 elements and 
lead are overestimated and would be better reproduced with a 
lower Mpmz or M\. 


The difficulty in reproducing the abundance ratios suggests 
that our simple parameterisation of the PMZ possibly ignores 
some important aspects of the physics involved in the problem. 
For example, the mass of our PMZ stays constant during the evo¬ 
lution of the AGB. In the neutron-capture process the 5-element 
peaks are filled progressively in time: in the first few thermal 
pulses the star produces mostly light-5 elements while in later 
thermal pulses the distribution is weighed more towards neutron- 
rich nuclei. A time-dependent Mpmz could modify the final light- 
s to heavy-.v ratios. For example, a PMZ that increases in time 
may be able to produce larger [hs/ls] and [Pb/hs] ratios. Values 
of [hs/ls] above 1 dex are observed in most CEMP-s/r stars and 
are currently impossible to reproduce in our models, as we will 
discuss in Paper II. On the other hand, a small PMZ in the last 
few pulses would probably enhance the light-5 peak compared 
to the heavy-5 peak and lead. 


The origin of the /'-element enrichment in CEMP-5 stars is an 
open issue. Several explanations have been proposed and some 
of these suggest that the r- and 5-enrichments in CEMP- 5 /r stars 
are independent (see Jonsell et al. 2006 1 . However, this interpre¬ 
tation does not explain the correlation observed in CEMP-5 stars 
between the abundances of barium and europium. 


6.4. The effect of non-convective mixing in dwarf CEMP stars 


Five stars in our sample have high surface gravities and therefore 
have not yet experienced the first dredge-up. The mixing mech¬ 
anisms that can modify the surface abundances in these stars are 
not well understood and many counteracting effects potentially 
play a role: thermohaline mixing, gravitational settling, radiative 
levitation and rotation. In our models we simulate two opposite 
situations: our default option is that the accreted material is di¬ 
luted throughout the entire star. Alternatively, the mixing is inef¬ 
ficient and the material remains on the surface until mixed in by 
convection. The choice of inefficient mixing improves the fit of 
two stars, CS22956-028 and CS29497-030. However, the fit of 
star CS22956-028 (Fig. |A.3| > is based on only five elements one 
of which, barium, is overestimated by 2 dex. An alternative in¬ 
terpretation of the observed abundances is that this star may have 
been polluted by a rotating massive star, which may produce up 
to [Sr/Ba] as 2 and [Pb/Sr] < —1 as suggested by the models of 
25 M q stars discussed by Frischknecht et al. ( ]20l 2[ ). The abun¬ 
dances of other elements, particularly nitrogen and lead, are nec¬ 
essary to test this hypothesis. On the other hand, CS29497-030 
is a CEMP- 5 /r star and the^ 2 of the fit is dominated by the large 
abundance of seven elements mainly produced by the /--process. 
A fit that includes efficient thermohaline mixing reproduces the 
5-process elements much better. In summary, our results do not 
provide strong constraints on the efficiency of non-convective 
mixing mechanisms in metal-poor dwarf stars. 
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6.5. The effect of fixed model metallicity 


The fact that we use a model tailored for [Fe/H] = -2.24 to in¬ 
vestigate stars of different iron abundances is an issue for many 
stars in our sample. By adopting a fixed metallicity for all our 
systems we implicitly ignore the effects due to the different 
[Fe/H]. This approximation likely introduces increasingly big¬ 
ger errors the larger the difference in [Fe/H]. Qualitatively, we 
expect the choice of a lower metallicity to result in two effects: 
an increase in the abundances of carbon and all the elements 
produced by AGB nucleosynthesis, because the iron abundance 
is smaller, and relatively larger abundances of neutron-rich iso¬ 
topes, because of the higher ratio of neutrons to iron seeds. If we 
vary the metallicity in our models, all other options being equal, 
we reproduce the first of the two effects but not the second, be¬ 
cause the neutron-to-seed ratio, and consequently the amount of 
i-process elements that are produced in our modelled intershell 
region, is determined by the mass of the star and by the mass of 
the PMZ at [Fe/H] = -2.24. Hence, we need to assume that the 
abundances of 5-isotopes in the intershell region do not change 
for metallicity variations. It is arguable which of the two options 
is the best: to adopt the same metallicity for all observed systems 
or to vary it in proportion to the observed [Fe/H]. In this study 
we preferred the first option because we are more confident of 
the results of our code at the fixed metallicity Z = 10 4 , but it is 
very important to keep in mind the limitations that derive from 
this choice. 


The most extreme example is CEMP-.v/r star CS29497-034 
that exhibits iron abundance [Fe/H] = -2.96 which corresponds 
to a metallicity five times lower than our default assumption. 
However, if we adopt Z = 10 4 we find a model that repro¬ 
duces the abun dances of all elements except europium (solid 
line in Fig. A.4 1 , while if we adopt Z = 2 x ICE 5 our best model 
overestimates the abundances of carbon, sodium and magnesium 
(dashed line in Fig. A.4 1 and does not improve the fit to the abun¬ 
dance of europium. This result suggests that the intershell com¬ 
position of AGB stars at [Fe/H] » -3 may be very different. 

Another example is HE0024-2523, that exhibits a distri¬ 
bution of 5-elements weighted towards neutron-rich elements 
and none of our model reconciles the large abundance of 
lead, [Pb/Fe] = 3.2, with the low enhancement of strontium, 
[Sr/Fe] = 0.5 (Fig. A.9 1 . This discrepancy can be qualitatively 
explained by the difference in the observed and modelled metal¬ 
licity (Z « 4 x 1CT 5 and Z = 10 4 , respectively). However, star 
CS22942-019 has almost the same iron abundance ([Fe/H] = 
-2.69) but we find a good match between our best model and 
the observations. 


7. Conclusions 

This study shows that the requirement to reproduce at the same 
time the chemical abundances and orbital periods observed in 
our sample of 15 CEMP-5 binary stars put strong constraints on 
the adopted binary evolution model. It is generally necessary that 
the modelled binary systems lose efficiently angular momentum 
and transfer mass with high accretion efficiency. In particular, 
binary systems with orbital periods below a few thousand days 
need to transfer mass more efficiently than normally assumed 
in our models. In a forthcoming paper we will analyse a larger 
sample of CEMP-5 stars without measured orbital periods. The 
comparison between the modelled periods determined in the two 
studies will have implications for our model of the wind mass- 
transfer process. 


When the condition of strong mass transfer is fulfilled the 
discrepancies between synthetic and observed abundances arise 
from the model of AGB nucleosynthesis. In about half of the sys¬ 
tems the observed element-to-element ratios are not reproduced. 
In particular, to match the large enhancements of the heavy-5 el¬ 
ements in CEMP-.v/r stars our models produce an excess of car¬ 
bon, sodium, magnesium and light-5 elements. This discrepancy 
and the fact that the abundances of the /--process elements are 
mostly underestimated suggest that in our model higher densi¬ 
ties of free neutrons should be produced in some circumstances, 
because a larger neutron-to-iron ratio would favour neutron-rich 
elements. 
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Appendix A: Best-fitting models to the observed 
CEMP-.S stars in our sample 

In this section we show the best-fitting models to the observed 
abundances of all CEMP-s stars in our sample. 
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atomic number 


A.l 


for CS22956-028. The best- 


Fig.A.3. Solid line and points: as Fig. 
fitting model assumes that non-convective mixing processes, such as 
thermohaline mixing, are inefficient and therefore the accreted material 
remains on the surface of the star. Dashed line: alternative model with 
efficient thermohaline mixing (x 2 = 51.8). 


atomic number 

Fig. A.l. Points with error bars: as Ftg.j^for BD+04 2466. Solid line: 
best-fitting model computed with model set A (Table[2]l. The combined 
abundance of carbon and nitrogen is shown in the right panel. 
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Fig. A.2. Star CS22948 - 027. Solid line and points: as Fig. 
the best model (with model set B) the secondary star accretes AM acc = 
0.27 Mq. With model set A (dashed line) the secondary star accretes 
AM acc = 0.12 M 0 and therefore the material is more strongly diluted. 
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Fig. A.4. Points with error bars: as in Fig. |A.l| for CS29497-034. Solid 
line: best-fitting model computed with set B adopting the default metal- 
licity Z = lO" 4 .^ i s dominated by the discrepancy in the abundance 
of europium (about 1 dex). If we exclude europium from the fit we find 
Tmi n = 15.2 and X^Jv = 1.7. Dashed line: alternative model computed 
with set B and reduced metallicity, Z = 2x 10~ 5 , that corresponds to the 
observed iron abundance [Fe/H] = -2.96. 
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Fig. A.5. As Fig. 
mens, such as N, Na, 
initial primary mass. 
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for CS29509-027. The abundance of other ele- 
Mg and Pb. is necessary to better constrain the 
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Fig. A.6. As Fig. A. 1 for HD198269. The best-fitting model (solid line ) 


computed adopting model set B. The abundances determined by Van 
|Eck et al7| < |2003|) and[Vanture (1992b) are represented as plus signs and 
filled circles, respectively. 
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Fig. A.7. As Fig. A. 1 for HD201626. The best-fitting model to the ob¬ 


served abundances is found with model set B. 
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Fig. A.8. As Fig. A.l for CEMP-j/r star HD224959. The abundances 


determined by [Van Eck et al.| ( |2003| and Vanture l |1992b| l are repre¬ 
sented as plus signs and filled circles, respectively. 
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for HE0024-2523. Solid 


Fig. A.9. Points with error bars: as Fig. 
line: Best-fitting model found with model set B and assuming a low effi¬ 
ciency for common-envelope ejection, q-ce = 0.03, which is required to 
reproduce the observed period of P orb = 3.14 days in a binary scenario. 
Dahsed line: alternative scenario in which the HE0024-2523 was ini¬ 
tially part of a hierarchical triple system. An intermediate-mass primary 
star was in a wide orbit around a close binary. During its AGB phase the 
primary star pollutes the inner binary. Subsequently, the secondary star 
overfills its Roche lobe, the system enters a common envelope the ejec¬ 
tion of which (modelled with or C E = TO) requires the orbit to shrink to 
the observed period. 
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Fig.A.10. As Fig. 
found with model set B. 
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for star HE0507 - 1430. The best-fit model is 
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Fig.A.ll. As in Fig. |A.l 
found with model set B. 
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for star LP625 - 44. The best-fit model is 
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Appendix B: Two-dimensional confidence intervals 

Fig. |B.l| t shows the two-dimensional confidence regions deter¬ 
mined for the initial orbital period, P„ and primary mass, M u, 
of model star CS29497-034. The one-dimensional confidence 
intervals of these two parameters are shown in Figs. B.l 3 and 
[BTj:, respectively, in which the same symbols as in Fig. 5] are 


used. Fig. B.2 is the same as Fig. B.l for star HD201626. To 


compute the two-dimensional confidence intervals of two input 
parameters p and p' we follow the same procedure as described 
in Sect. 4.3 but we fix a pair of values (p, p') and we let the other 
two parameters vary. In two dimensions the confidence levels of 
68.3%, 95.4% and 99.7% correspond respectively to the thresh¬ 
olds Ax 2 = 2.30, 6.17, 11.8 (Press et al. 1989). These thresh¬ 


olds are shown in Figs. B. 1 1 and B.2 1 as solid, dashed and dotted 
lines, respectively. 

The two-dimensional confidence intervals in Fig. |B.l| i indi¬ 
cate that there is correlation between the initial orbital period and 
primary mass of the model systems. The more massive the pri¬ 
mary star i s, the longer needs to be the initial period to minimise 
X 2 - In Fig. B.2i the correlation between the two parameters in 
the models of star HD201626 is even more clear. This correla¬ 
tion indicates that although the confidence range of M\ , is rather 
large, as shown in Table [4] for each primary mass in our grid 
there is essentially only one orbital period at which the model 
reproduces the observations. 





Fig. B.l. Confidence intervals of the input parameters for model 
star CS29497-034. Panels b and c show the one-dimensional 
confidence intervals of P t and M\\, respectively. Long-dashed 
lines indicate the thresholds Ax 2 = 1, 4, 9. Panel a shows the 
two-dimensional confidence intervals. The thresholds Ax 2 = 
2.30,6.17, 11.8 are represented as solid, dashed and dotted 
lines, respectively. The blue plus sign indicates the best model. 




Fig. B.2. Confidence intervals of the input parameters for model 


star HD201626. The symbols are as in Fig. B.l 





































